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Abstract: We consider p- variations in some membrane potential data -viewed as a function of the step size in 
case where p is fixed, or viewed as a function of p in case where the step size is fixed- and compare their shape 
with results in Jacod and Ait-Sahalia [1] which do hold for general semimartingales. We obtain the following 
conclusion: in non- or very rarely-spiking cases the membrane potential behaves as a semimartingale, in some 
cases as a semimartingale with jumps. Once the neuron is spiking, a semimartingale modelization is no longer 
adequate for the membrane potential between successive spikes, even if interspike intervals are relatively long. 
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We take a new look on two data sets recording the membrane potential in a pyramidal neuron (intra- 
celluar recording) which belongs to a cortical slice observed in vitro (representing an active network) 
under different experimental conditions. The neuron under observation receives synaptic input from a 
large number of other neurons in the slice. Stimulating the slice -and thus the networking properties 
of all neurons belonging to the slice- by a potassium bath, W. Kilb (Institute of Physiology, University 



of Mainz) recorded 'Zelle 3' in 2004 and '17Sept08_023' in 2008. The data are shown in figures [20 
and [21] In 'Zelle 3', 10 different concentrations of potassium correspond to 10 different data sets 
(called 'levels' below, obtained under 3, 4, 5, 6, 7, 8, 9, 10, 12, 15 mM of K) observed over 60 seconds 
each. In '17Sept08_023', one potassium level (5 mM of K) was kept constant over a much longer time 
interval; in the present note, for ease of comparison with 'Zelle 3', we use only the first 60 seconds of 
observation from this data set. 

In terms of a diffusion process modelization, 'Zelle 3' has been considered in [4], '17Sept08_023' in [6] 
(section 5.3 there, using the estimation method of [4]). In [3] and [6], assuming that the membrane 
potential between successive spikes (more precisely: sufficiently away from the spikes) can be modelled 
as a time homogeneous diffusion process, nonparametric estimates for diffusion coefficient and drift 
made appear a linear mean-reverting drift combined with either a constant or a linear or a 'bowl- 
shaped' diffusion coefficient: these cases correspond to Ornstein-Uhlenbeck (OU) type, Cox-Ingersoll- 



Zelle 3 K 15.asc, 60 sec, 2-variations depending on M 
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Figure 1: Membrane potential 'Zelle 3' level 10 (15 mM of K): plotting truncated 2-variations M -> Vr(2, A, M) 
as defined by ^ for 1 < M < 240, in increasing order for T = 4, 8, 10, 16; no further changes above V = 10. 

Ross (CIR) type or -in the language of |3|- Pearson (P) type diffusions. Analyzing the data in the 
same way in smaller time windows, we can assert that the assumption of time homogeneity seems well 
satisfied in '17Sept08_023' (here W. Kilb had used a new type of electronic stabilization device), and 
reasonably well satisfied in several potassium levels of 'Zelle 3'; obvious exceptions are the 'low' levels 
1, 2, 4 where strong time inhomogeneities appear (discussed for level 1 in [3], section 4.6). 

In this note, we consider for p=2 or p=4 fixed p-variations in the membrane potential data as a 
function of the step size, i.e. the length of the time intervals over which variations are calculated, and 
compare these to simulated diffusion equivalents whose drift and diffusion coefficients are as estimated 
in [3] or [6] . Relying on recent results of Ait-Sahalia and Jacod [1] , we then ask the question to which 
extent a (continuous) semimartingale model is in fact adequate for the membrane potential between 
successive spikes. In our data, a surprising difference appears -and in particular in the same neuron 
'Zelle 3'- between spiking and non-spiking regimes. We then fix the step size and consider p- variations 
in our data as a function of p: again the same striking difference between spiking and non-spiking 
regimes arises. When spiking is sufficiently frequent (in our data, in levels 9 and 10 of 'Zelle 3', and 
in '17Sept08_023'), pictures of power variations arise which do not agree with what we should see in a 
semimartingale, continuous or not, according to [1J. This is not simply an effect of noisy observation 
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Zelle_3_K_9.asc. whole trajectory, 2-variations depending on M 
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Figure 2: Membrane potential Zelle 3 level 7 (9 mM of K): plotting truncated 2-variations M — !• V r (2, A, M) 
as denned by in increasing order for T 6 {1, 2, 4, 8, 10, 16, 32}; no further changes above F = 8. 



Zelle_3_K_9-simul_aequiv, whole trajectory, 2-variations depending on M 
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Figure 3: Simulated P type diffusion equivalent (jl], section 4.4) for Zelle 3 level 7: 2-variations plotted in 
analogy to figure [2j 
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Zelle_3_K_9.asc, whole trajectory, 4-variations depending on M 
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Figure 4: Membrane potential Zelle 3 level 7 (9 mM of K): plotting truncated 4-variations M — !• Vr(4, A,M) 
as denned by in increasing order for T 6 {1, 2, 4, 8, 10, 16, 32}; no further changes above F = 8. 



Zelle_3_K_9-simul_aequiv, whole trajectory, 4-variations depending on M 
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Figure 5: Simulated P type diffusion equivalent for 'Zelle 3' level 7 ([?], section 4.4): 4-variations plotted in 
analogy to figure |4j 
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(equally observable in all levels of 'Zelle 3', but almost absent in '17Sept08_023' where a different type 
of electrode had been used) but concerns the shape of the curve of power variations as a function of 
the step size when p is fixed, or as a function of p when the step size is fixed. In sharp contrast to 
this, in non- or rarely spiking regimes (the spikeless levels 1-7 of 'Zelle 3' (3-9 mM of K), and level 
8 (10 mM of K) where one single spike is emitted during the overall observation time of 60 seconds) 
the pictures of power variations agree very well with what is to be expected for a semimartingale -up 
to secondary effects like noisy observation or feedback effects in the slice- and with what can be seen 
in simulated diffusions or jump diffusions. 

Our conclusion is that a semimartingale model seems adequate for neurons in non-spiking or 
rarely-spiking regimes, whereas something essentially different -not well captured by semimartingale 
modelization- seems to prevail in spiking regimes. 

This note is organized as follows. Section 1 considers for fixed p (p = 2 or p = 4) p-variations 
as a function of the step size over which we calculate the increments. Subsection 1.1 explains the 
truncated power variations which we use in this note, subsection 1.2 considers the spikeless or very 
rarely spiking levels of 'Zelle 3', subsection 1.3 the spiking levels of 'Zelle 3' and the frequently spiking 
neuron T7Sept08_023'. Section 2 considers p-variations as a function of p for fixed step size, with 
an analogous program. I would like to stress that this note is a 'not really mathematical' paper (no 
theorem, no rigorous proofs, some merely 'plausible' approximations): its aim is to analyze a set of 
neuronal data in the light of theorems in Ait-Sahalia and Jacod pQ which do hold for very general 
semimartingales, and to show that some essential difference exists between spiking and non-spiking 
regimes (in the same neuron) in view of semimartingale modelization. I would like to thank H. 
Luhmann and W. Kilb for the data, and J. Jacod for some longer discussions on this problem. 



1. Fixing the power p and varying the size of the increments 

The structure of the data is as follows. The different experiments in 'Zelle 3' record membrane 
potentials at times tj := iA, < i < 100001, with A = 6 • 10~ 4 [sec], thus with total observation time 
T = 60 [sec]. The data '17Sept08_023' are considered only in restriction to the first 60 seconds of 
observation: here the time grid is U = iA with A = 2 • 10~ 4 [sec], < i < 300001. The measurement 
is in millivolt [mV], formally with three decimals, but with an information that the third decimal 
is not reliable at all. When spikes are present in the data set, we remove time neighbourhoods 
(r — 0.12, r + 0.18) [sec] centred at the spike times r from the data, in order to exclude any influence of 
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Figure 6: Replacing dW t in the simulated diffusion equivalent to 'Zelle 3' level 5 ([3], section 4.1) by dW t +edS" 
for a = 1.75 and e — 0.1: truncated 4-variations M —> Vr(4, A,M) as denned in ^ for 1 < M < 10, in 
increasing order for T E {1, 2, 4, 8, 10, 16, 32, 64, 128, 256}; no further changes above T = 64. 

the typical shape of the spike. We calculate increments, p-variations, . . . over spikeless time segments, 
and then add up corresponding terms coming from different segments. The estimators which we use for 
drift and diffusion coefficient in a discretely observed diffusion are those of [1], and are not explained 
here. Here we explain the way we calculate a p- variation over a time segment, and over a collection of 
time segments, in view of application to the membrane potential in a neuron which can emit spikes. 
In the present section, we concentrate on fixed power p and vary the step size, i.e. the length of the 
time intervals on which increments are evaluated. 

1.1. Truncated p- variations for neuronal data 

First, for varying choices of a truncation factor < T < oo and for multiples M of the step size A 
prescribed by the data, for p > 2 fixed, we define 

, h-M 

(1) Vt , tl ,v(p, A, M) := — \ X (i+M)A- X iA \ P \\ X(i+M)A -X iA \ <3VAMT} 

i=io 

with respect to one spikeless segment [to>^i] = [*oA, iiA] (we define a spikeless segment as a maximal 
interval between r r „i + 0.18 and r r — 0.12, avoiding neighbourhoods of the successive spike times r r _i, 
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defined above). Second, based on ([!]), we define 

(2) V T (p,A,M) := V t0il>tlAjT (pA,M) + ... + V t0ii , tlA , r ( P ,A,M) 

for the whole membrane potential trajectory up to time T = 60 [sec], where [to,l>*l,l]> ••• > [to,i,h t e] 
denotes the collection of spikeless segments (including an initial [0, t\ — 0.12] before the first spike and 
a final [rg-i + 0.18, T] after the last spike). For spikeless membrane potentials or for the simulated 
diffusion equivalents, the full interval [0, T] is the unique segment. 

With truncation factor T increasing to oo in ([2]), we will finally capture all jumps of a semimartingale 
trajectory up to time T, or all increments from spikeless segments in a membrane potential data set, 
and will arrive for V tending to oo at 

it—M 

M 



h-M 

(3) V khtl (p,A,M) := V t(htliOQ (p,A,M) = - ^ \X 



(i+M)A — -XiA I 



1=10 



for a single spikeless segment [to,*i]) arid at 

(4) V(p,A,M) := V t0tUtlil (p,A,M) + ... + Vt 0tt , tltt (p,A,M). 

for the whole membrane potential trajectory. For to < s o < t±, consider MA-step p- variations on 
[so^i] as defined in Ait-Sahalia and Jacod ([Tj, formula (9)): 

L MA J 

Bs ,ti(p,A,M) := \X So +kMA - ^s +(fc-l)A/A| P 

fc=l 

and note that the right hand side of ^ equals 

. M-l 

(5) m £%, +J - )A>tl (p,A,M). 

3=0 

Averaging over j = 0, 1, ... , M—l in ^ allows to make use of all MA-step increments available in 
the time window [ig,ti]. Heuristically, for MA sufficiently small, all summands in §5§ should be very 
close to Bt lh t! (p, A, M), thus we will make the following approximation during the present note : 

. M-l 

(6) ^ X)B( <0+i )A,t 1 (p,A,M) « B Uhtl (p,A,M) . 

Since the time resolution A of the data set cannot be modified by the statistician, asymptotic results 
as given in Jacod and Ait-Sahalia (p], (11)-(13)) for time-step tending to have to be mimicked 
through variation of multiples M of A. Assuming that MA is sufficiently small for the M which we 



Zelle_3_K_12.asc. whole trajectory, 2-variations depending on M 
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Figure 7: Membrane potential 'Zelle 3' level 9 (12 mM of K): plotting the 2-variations M Vr(4,A,M) as 
denned by j2j, in increasing order for V £ {1,2, 4, 8, 10, 16, 32}; no further changes above T = 8. 



Zelle_3_K_12-simul_aequiv, whole trajectory, 2-variations depending on M 



40 

range(Mwerte) 



Figure 8: Simulated CIR type diffusion equivalent for 'Zelle 3' level 9 ([3], section 3.3): 2-variations plotted in 
analogy to figure [7J 
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Zelle_3_K_12.asc, whole trajectory, 4-variations depending on M 
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Figure 9: Membrane potential 'Zelle 3' level 9 (12 mM of K): plotting the 4-variations M -t Vr(4, A, M) as 
denned by in increasing order for V £ {1, 2, 4, 8, 10, 16, 32}; no further changes above T = 8. For T > 8, we 
calculate values 8.17 for M = 3, 5.42 for M = 2, 3.20 for M = 1. 



Zelle_3_K_12-simul_aequiv, whole trajectory, 4-variations depending on M 




20 40 60 

range(Mwerte) 



Figure 10: Simulated CIR type diffusion equivalent for 'Zelle 3' level 9 ([3], section 3.3): 4-variations plotted 
in analogy to figure [9] 
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consider, and that the data in fact do correspond to a discretely observed semimartingale £ = (£ s ) s >o, 
we may read (pQ, (ll)+(7)+(10)) as 

(7) B t0>tl (4,A,M) « B t0>tl (4) := ^ |A£ S | 4 as M = . . . , 3, 2, 1 gets small 

(8) B t0>tl (2,A,M) « A t0)tl (2) + B t0)tl (2) := a 2 (£ s )ds+ |A£ S | 2 as M gets small 

(9) for £ continuous : Bt 0)tl (4, A, M) « M • 5to, tl (4, A, 1) as M gets small . 

By ([T]), 4- variations stabilizing at a strictly positive 'limit' when M gets small indicate presence of 
jumps in the semimartingale £. For £ continuous, 4-variations should be linear in M as long as M 
is small, as a consequence of (g). On every segment, by ([I], theorem 1), this is a dichotomy which 
represents a test for presence of jumps in a semimartingale £ = (£t)t>o recorded at time resolution A. 
Putting together the segments as in Q above, we rephrase the test of [I] in the following form: as 
M = ... 3, 2, 1 gets small, 

(10) for £ with jumps : V{A, A, M) stabilizes at a strictly positive 'limit' ; 

(11) for £ continuous : M — > V(4, A,M) is linear . 



The results of Ait-Sahalia and Jacod [T] being asymptotic results for shrinking time grids on which the 
process is observed, reformulations such as ([7])+(|9]) or (10)+(11) of this test hinge on the assumption 
that MA be 'sufficiently small' for the M which we wish to consider. In the data, we can not modify 
the time resolution A. It may well happen that considering MA for 1 < M < 5 (say), we are not 
yet 'sufficiently small' in the sense of [I]. As an example, replace in the OU diffusion equivalent to 
'Zelle 3' level 3 (as in @J, section 4.1) the driving dW t by dW t + edS? for small e where S a = (S?) t > 
is a symmetric stable process with index a £ (0,2). Simulating increments of S a using Chambers, 
Mallows and Stuck [2], with A the time resolution of 'Zelle 3', the test ([7])+([9]) will be unable to detect 
presence of jumps in the simulated jump diffusion for a very close to 2, whereas in case a = 1.75, the 
jumps are detected (see figure pi) by inspection of 4-variations for M < 5. 



1.2. Application to the non- or very rarely spiking levels of 'Zelle 3' 

If we admit heuristics ([6]), then figures [2]+|4] (2-variations and 4-variations for level 7 of 'Zelle 3', no 
spikes) in comparison to figures [3]+[5] (2-variations and 4-variations for a simulated P type diffusion 
equivalent, with drift and diffusion coefficient as estimated in [4], section 4.4) show that level 7 of 
'Zelle 3' exhibits the typical features of a continuous semimartingale, up to some strong deformation 
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of the initial part of the 2- variation as a function of M, visible for M- values up to ~ 10. We interpret 
this deformation as noise contaminating the observation, generated by the electrode measuring the 
membrane potential. This is supported by the observation that in all levels 1, . . . , 10 of 'Zelle 3' the 
2-variations attain a value close to ~ 300 for M = 1 (see figures [T]+[7]+[2]) , whereas in the recording of 
T7Sept08_023' a different type of electrode was used which does not produce the same phenomenon 



(figure 11 ). 



Qualitative agreement (disregarding the effect of noise for M < 10) between 2- and 4-variations for 
the neuronal data and 2- and 4-variations for their simulated diffusion equivalents is observed in all 
spikeless levels 1, . . . , 7 of 'Zelle 3', and also in level 8 where one single spike is generated in the 60 
seconds of observation. This agreement is not always as perfect as in figures [2]+[3] and [4]-|-[5] (presenting 
level 7 of 'Zelle 3' in the figures, we did chose the level where the best fit occurred), but the qualitative 
features (2-variations flat in M up to the initial effect of noise, 4-variations linear in M) agree well 
between data and simulated diffusion equivalent. Hence in the spikeless or rarely spiking levels of 
'Zelle 3', the membrane potential (away from the isolated spike in case of level 8) can be be viewed as 
a semimartingale. 

Is the semimartingale continuous, or does it have jumps? The 4-variations for levels 1-8 look at 
first glance very much like being linear in M for small M-values. In some of these levels however, a 
closer look to small M values might suggest presence of jumps. As an example, comparing figures 
|4]+[5] (level 7) for M = 3,2, 1, there is a difference in the behaviour for small M, and in figure [4] we 
may see convergence to some strictly positive 'limit' as M gets small. Thus there might be jumps 
in the membrane potential data 'Zelle 3' level 7. The same effect is visible e.g. in levels 1 (3 mM of 
K), 3 (5 mM of K), 4 (6 mM of K). However, since in all pictures of 2-variations we saw additional 
variation for small values of M, with the interpretation of noise of the measuring electrode, this noise 
might similiarly affect the 4-variations for small values of M. Hence, with the methods of section 1, we 
cannot decide whether or not the non- or rarely spiking levels of 'Zelle 3' should be viewed as continuous 
semimartingales or as semimartingales with jumps. We will be able to answer this question with a 
different method in section 2. 

1.3. Application to the spiking levels of 'Zelle 3', and to the neuron '17Sept08 023' 

In the spiking levels 9 and 10 of 'Zelle 3' (18 spikes in level 9, 8 spikes in level 10, over a total of 
60 seconds of observation time) and in the rapidly spiking neuron T7Sept08_023' (~ 50 spikes over 

in 



the first 60 seconds), pictures very different from those discussed above arise (see figures 11 + 13 
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17Sept08_023.asc, on [ 0, 60 ], 2-variations depending on M 




M ( using increments over M * 2e-04 sec, with 2e-04 the time resolution in data set) 

Figure 11: Membrane potential '17Sept08_023' (5 mM of K): plotting 2-variations M V r {2, A, M) as defined 
in ([2]), with truncation constant T, in increasing order for T 6 {1, 2, 4, 8, 10, 16, 32, 64, 128}; no further changes 
above T = 64. 

17Sept08_023-simul_aequiv, 60 sec, 2-variations depending on M 




M ( using increments over M * 2e-04 sec, with 2e-04 the time resolution in data set) 



Figure 12: Simulated CIR type diffusion equivalent ([BJ, section 5.3) for 17Sept08_023 : 2-variations calculated 
in analogy to figure 11 
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17Sept08_023.asc, on [ 0, 60 ], 4-variations depending on M 
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M ( using increments over M * 2e-04 sec, with 2e-04 the time resolution in data set) 



Figure 13: Membrane potential '19Sept08_023' (5 mM of K): plotting 4-variations M -> V T (A, A, M) as denned 
in ([2]) with truncation constant T, in increasing order for T 6 {1, 2, 4, 8, 10, 16, 32, 64, 128}; no changes for T > 64. 



1 7Sept08_023-simul_aequiv, 60 sec, 4-variations depending on M 
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M ( using increments over M * 2e-04 sec, with 2e-04 the time resolution in data set) 



Figure 14: Simulated CIR type diffusion equivalent ([5J, section 5.3) for '17Sept08_023': 4-variations calculated 



in analogy to figure 13 
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comparison to |T2|+ 14 for '17Sept08_023', see figures [7|+(9] in comparison to [8]+ 10 for 'Zelle 3' level 9, 



see also figure [l]for 'Zelle 3' level 10). There is no longer a qualitative coincidence between the shape 
of the p- variations as a function of M in the data, and the shape of the p- variations as a function of 
M in the simulated diffusion equivalent. 

For the 4-variations of '17Sept08_023' (where thanks to A = 2 • 10~ 4 [sec] we get nearer to than 
in the data 'Zelle 3') we calculate the values 39.13 at M = 3, 12.50 at M = 2, 1.73 at M = 1 for 



sufficiently large truncation factor T, see the detail in figure 15 Under p = 4, these values rule out at 
the same time 



the possibility of a strictly positive 'limit' in (10) when M gets small 



the possibility of a linear dependence on M in (11) when M gets small 



and thus -according to the dichotomy in (|TJ, theorem 1)- rule out the possibility that the data 
T7Sept08_023' represent a discretely observed semimartingale £ = (£t)t>o- With respect to this result, 
problems such as time inhomogeneity, obviously present in at least a part of our data, or presence of 
jumps as raised in [JJ become irrelevant. Curves of similiar shape are obtained for the 4-variations in 
'Zelle 3' level 9 (figure [9j and level 10. The situation for small values of M is less clear in the data 
'Zelle 3' than in T7Sept08_023' (the 4-variations in 'Zelle 3' level 9 take the values 8.17 at M = 3, 
5.42 at M = 2, 3.20 at M = 1, cf. figure |9| and may be affected by the noise which appears clearly 
in the 2- variations for small M). In all three data sets, the shape of the 2-variations as function of 
the step size (figures [7J [I 11), with a remarkable maximum at M ~ 80 for 'Zelle 3', at M ~ 40 for 



'17Sept08_023', is incompatible with a discretely observed semimartingale, cf. ([8]). Thus we conclude 
that the membrane potential in the spiking neuron -between successive spikes, staying sufficiently 
away from the spikes- is not adequately modelled by a semimartingale. 

We remark that despite this fact, nonpar ametric estimates for drift and diffusion coefficient -within a 
semimartingale setting- may produce seemingly satisfactory results. As an example, for 'Zelle 3' level 
10, one obtains a convincing fit between occupation time calculated from the data and the invariant 
Gamma type law of the CIR type diffusion equivalent for 'Zelle 3' level 10 (with estimated drift and 
diffusion coefficient according to [3], section 3.2). But also here, the hypothesis of a diffusion process 
became questionable when in case of 'Zelle 3' level 10 the estimates used in [4] were observed to depend 
much more on the chosen multiple M of the step size A -entering the definition of the kernel estimator 
in [4] - than was claimed in [3j. This observation represents a surprising contrast to what has been 
checked for the non-spiking levels 3, 6, 7 of 'Zelle 3' in (pi], figure 10). 
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17Sept08_023.asc, on [ 0, 60 ], 4-variations depending on M 




M ( using increments over M * 2e-04 sec, with 2e-04 the time resolution in data set) 



Figure 15: Zooming into figure 13 for small values of M: for T > 64, we obtain the values 39.13 for M = 3, 
12.50 for M = 2, 1.73 for M = 1. 



1.4. An additional remark 

In all levels of 'Zelle 3', the 2-variations M — > Vr(2, A, M) present periodic deformations, for fixed 
value of the truncation factor T; to less extent, this is visible also in the 4-variations M — > Vr(4, A, M) 
(see figures [2]+|4j (7[f[9j [T]) . Independently of the level and of the value of T, the deformations are most 
visible near M ~ 32 and M ~ 64 (in figure [T] continued by ~ 96, 128, . . .), and go attenuating as 
M gets larger. This might indicate that the neuronal network in the slice to which 'Zelle 3' belongs 
possesses loops or circuits, and thus produces feedback at a fixed periodicity. 

2. Fixing the size of the increments and varying the power p 



In section 1, we have considered p- variations for fixed pas a function of the step size. We continue this 
discussion, but now with p- variations considered as a function of p for fixed step size. Assuming that the 
spikeless segments of the membrane potential do correspond to a discretely observed semimartingale 
with nonvanishing continuous local martingale part (an unproblematic additional assumption given 
the shape of the power variations for small truncation factor T appearing in figures [2j |4j [7j [9j [TTJ [13J 
15 above) and that both M' £ {M,2M} lead to sufficiently small values of M'A, we may read (PQ, 
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Zelle_3_K_3.asc, 60 sec, test for jumps through p-variations functional in p 




trunc. factor = 4 



trunc. factor - 886 



12 3 4 

p ( using increments over 2M / M times the time resolution 6e-04 [sec] in the data set, M = 20 ) 



Figure 16: Membrane potential 'Zelle 3' level 1 (3 mM of K, no spikes). Logarithm of ratios y^f^TMy as a 
function of p, for V G {1,4,16,64,256}. Increasing values of the truncation factor correspond to more 'solid' 
red curves; no changes occur above T — 16. It is seen that the data 'Zelle 3' level 1 are well compatible with a 
continuous semimartingale observed at discrete times «A, < i < 100001. 



(ll)+(7)+(10)) on a segment [to>*i] as follows: 



(12) for £ with jumps : B t0jtl (p, A, 2M) 



(13) for £ continuous : B t0ttl (p, A, 2M) 



B t0}tl (p,A,M) for2<p<oo 
2 E 2- 1 B t(htl (p,A,M) for0<p<2 
25- 1 B toM (p, A, M) for < p < oo 



Again we accept the heuristics of section 1.1, in particular the approximation ([6]). We extend the 
heuristics by assuming that if spikeless segments of the membrane potential do correspond to a semi- 
martingale which has jumps, then jumps will occur on every segment [to, t±] under consideration (this 
is unproblematic e.g. if the Levy measure of the jump part of £ has infinite total mass independently 
of time). Then we can rephrase the test for jumps in Ait-Sahalia and Jacod (pQ, theorem 1) for fixed 
M and varying p as follows: 



(14) for £ with jumps : p 



, V(P,A,2M) . . 
log — — — - — — — - is approximately < 



(15) for £ continuous : p — > log 



V( P ,A,M) 

V(p,A,2M) 
V(p,A,M) 



constant = on [2,oo) 
linear in p on (0, 2) 

is approximately linear in p on (0, oo) . 
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Figure 17: Membrane potential 'Zelle 3' level 4 (6 mM of K, no spikes): logarithm of ratios a m) as a 
function of p, for V G {1,4,16,64,256}. Increasing values of the truncation factor correspond to more 'solid' 
red curves, no changes above V = 16. It is seen that the data 'Zelle 3' level 1 are well compatible with a 
semimartingale which has jumps, observed at discrete times zA, < i < 100001. 



The slope of the linear parts in (14)+(15) is deterministic, by (12)+(13). Thus, from the very begin 



ning, we know the shape which we expect to see when the membrane potential data between successive 
spikes do correspond to a discretely observed semimartingale (^t)t>0 : hi this case, the empirical object 

V(p, A, 2M) 



P 



log 



V(p,A,M) 



in (14)+(15) should be close to 



(16) 
(17) 



P 
P 



min | ( ^ — 1) log 2 , o| on < p < oo if £ has jumps , 
P 

(— — 1) log 2 on < p < oo if £ is continuous . 



Comparing the empirical object in (14)+(15) to the truncated line (16) expected for a semimartingale 



having jumps, or to the straight line (17) expected for a continuous semimartingale, we can decide 



whether jumps are present. Beyond this, since the above (14)+(15) represents a dichotomy on a very 



general class of semimartingales, see [T], we can decide whether or not our membrane potential data 
(away from the spikes) do correspond to a semimartingale. 



When visualizing the empirical object in (14)+(15), we will continue to make use of the truncation 
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17Sept08_023.asc, 60 sec, test for jumps through p-variations functional in p 




12 3 4 

p ( using increments over 2M / M times the time resolution 2e-04 [sec] in the data set, M = 20 ) 



Figure 18: Frequently spiking neuron '17Sept08_023' (5 mM of K, ps 50 spikes over the first 60 seconds): 
logarithm of ratios ^^(pA™) Pitted as a function of p, for T G {1,4,16,64,256}; the curves stabilize for 
r > 64. Increasing values of the truncation factor correspond to more 'solid' red curves. The dotted black lines 
show what is expected for a semimartingale, continuous or not, by [1]. For T tending to oo, the logarithmic 
ratios calculated from the data '17Sept08_023' turn out to be quite far away from a semimartingale hypothesis. 



factor r as in Q+Q, and calculate from our data log-ratios 

(18) ' ~> '° g V r(p ,A,M) ■ 

for varying values of V which are representative for asymptotics r — ¥ oo. 

Application to the data 'Zelle 3' and to '17Sept08_023' 



Figure 16 shows the lowest level of 'Zelle 3' (3 mM of K, spikeless): the picture corresponds very well 



to what we expect for a semimartingale which is continuous, as explained in (16)+(17)+(18) above. 



Figure 17 shows the level 4 of 'Zelle 3' (6 mM of K, spikeless): the shape of the curve corresponds very 
well to what we expect for a semimartingale which has jumps. Simulated diffusion equivalents (resp.: 
simulating a jump diffusion as in figure ml in relation to 'Zelle 3' level 5) produce pictures similiar 



to figure [16] (resp.: to figure 17). Moreover, all non-spiking levels 1-7 of 'Zelle 3', and even level 8 
with one isolated spike over 60 seconds of observation, lead to curves corresponding convincingly to a 
semimartingale hypothesis. Among these, exactly two -the levels 4 and 5 (level 5, not shown, looks 
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Figure 19: 'Zelle 3' level 10 (15 mM of K, 8 spikes over 60 seconds): logarithmic ratios am) plotted as 
a function of p, for T e {1,4, 16,64,256}; the curves stabilize for V > 64. Increasing values of the truncation 
factor correspond to more 'solid' red curves. The dotted black lines show what is expected for a semimartingale, 
continuous or not, by [I]. For T tending to oo, the logarithmic ratios calculated from the data T7Sept08_023' 
do not fit well with a semimartingale hypothesis. 



much like figure 17)- indicate the presence of jumps. Thus the method used in the present section, in 



contrast to the method used in section 1, is able to answer the problem of jumps raised in section 1. 

It turned out in section 1 that the membrane potential in the spiking levels 9 and 10 of 'Zelle 3' and in 
the frequently spiking neuron T7Sept08_023' was not adequately modelled by a semimartingale. The 



method of the present section reinforces this, see figures 18 + 19 The picture for level 9 of 'Zelle 3' (not 



shown) is similiar to what we show for level 10 in figure 19 The three curves do not correspond to what 
we expect for a semimartingale -continuous or not- by pp. Note that in levels 9 or 10 of 'Zelle 3', the 
time intervals between successive spikes are still relatively large and contain enough observations such 
that semimartingale methods -if the observed process were a semimartingale- should work successfully. 



On the basis of (16)+(17), we conclude as in section 1 that in the spiking levels of 'Zelle 3' and in 
the frequently spiking neuron T7Sept08_023', the membrane potential between successive spikes is not 
a discretely observed semimartingale. Note that our data 'Zelle 3' have been collected in the same 
neuron in the same cortical slice under different level of stimulation where stimulation by potassium 
activates the networking properties of all neurons in the slice. The membrane potential of the observed 
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neuron -sufficiently away from the spike times whenever there are spikes- behaves as a semimartingale 
as long as there are no spikes or at most extremely isolated ones, and loses this property once spikes 
occur frequently enough. This adresses in particular a serious question to some widely used neuronal 
models where interspike intervals are identified with level crossing times of semimartingales. 
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'Zelle_3': one pyramidal neuron, cortical slice, potassium (K) bath of Increasing concentration 




Figure 20: 'Zelle 3', all levels 1-10: membrane potential in the same pyramidal neuron under different exper- 
imental conditions. The neuron belongs to a cortical slice observed in vitro. The networking properties of all 
neurons in the slice are stimulated by a potassium bath (3, 4, 5, 6, 7, 8, 9, 10, 12, 15 mM of K). Spikes occur 
in levels 9 (18 spikes within 60 seconds of observation time) and level 10 (8 spikes), one isolated spike being 
observed in level 8. The time resolution is A = 6 • 10~ 3 [sec]. Data from H. Luhmann and W. Kilb, Institute 
of Physiology, University of Mainz. 
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Figure 21: Membrane potential in the frequently spiking neuron T7Sept08_023' (5 mM of K); we will use 
only the part of the data which corresponds to the first 60 seconds of observation. On this time interval, 
approximately 50 spikes occur. The time resolution is A = 2 • 10 -3 [sec]. Data from H. Luhmann and W. Kilb, 
Institute of Physiology, University of Mainz. 
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